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ABSTRACT 

Advances in theoretical ideas on how galaxies formed have not been strongly influenced by the advances 
in observations of what might be in the voids between the concentrations of ordinary optically selected 
galaxies. The theory and observations are maturing, and the search for a reconciliation offers a promising 
opportunity to improve our understanding of cosmic evolution. I comment on the development of this 
situation and present an update of a nearest neighbor measure of the void phenomenon that may be of 
use in evaluating theories of galaxy formation. 
Subject headings: cosmology: theory — galaxies: formation 



1. INTRODUCTION 

Voids between the regions occupied by normal textbook 
galaxies contain few observable galaxies of stars or gas 
clouds, and, apart from a tendency to greater gas content 
and star formation rates, objects observed in and near 
voids seem to be close to a fair sample of the cosmic mix. 
This striking effect will be termed the void phenomenon. A 
numerical definition, based on nearest neighbor statistics, 
is summarized in §3.4. 

The void phenomenon may be related to the preferences 
of early-type galaxies for high density regions and of galax- 
ies rich in gas and star formation for the edges of voids. 
This morphology-density correlation follows in a natural 
way from the biased galaxy formation picture in the Gaus- 
sian adiabatic cold dark matter (CDM) model. Perhaps 
the biased formation extends to the voids, where the mor- 
phological mix swings to favor dark galaxies. But if this 
were so the shift to the void mix would have to be close 
to discontinuous (§2.2), in remarkable contrast to the ob- 
served relatively slow variation of morphological mix with 
ambient density in regions occupied by visible galaxies. 
From continuity one might have thought the more likely 
picture is that gravity has emptied the voids of mass as 
well as galaxies. This does not happen in the CDM model, 
however. Simulations show, between the concentrations of 
large dark mass halos, clumps of mass that seem to be ca- 
pable of developing into void objects observable as clumps 
of stars or gas, contrary to what is observed. 

The apparent discrepancy has not attracted much atten- 
tion since the introduction of the biased galaxy formation 
picture fifteen years ago. This is partly because we do not 
have an established theory of how mass concentrations be- 
come observable (though we do have some guidance from 
the observations of galaxies at ambient conditions that do 
not seem to be very different from the voids; §4.1). A con- 
tributing factor is that the theoretical community has not 
settled on a standard numerical measure of the void phe- 
nomenon. A nearest neighbor statistic is commonly used 
in the observational community. This is discussed in §3, 
and §3.4 summarizes an update of the nearest neighbor 
statistic that might be applied to realizations of galaxy 
formation models. 

The main points of this paper are summarized in §5.1. 



Section 5.2 offers comments on the history of the void phe- 
nomenon as an example of the complex and occasionally 
weak interaction of theory and practice in a developing 
subject like cosmology. The discussion of the observational 
situation in §2 and §3 and of the theoretical situation in §4 
illustrates the unusually lengthy duration of modest inter- 
action between theoretical and observational ideas about 
voids, and the resulting opportunity now to learn some- 
thing new. 

2. VOIDS AND VOID OBJECTS 
2.1. Voids 



Rood (1981) gives an excellent picture of the state of ob- 
servational studies of galaxy clust ering tw o decades ago. 
He remarks that 20 years earlier Mayalj (1960, fig. 3) 
had found that in the direction of the Coma cluster far 
more galaxies are in the cluster than the foreground. Red- 
shift surveys to l arger angular distances fr om th e cluster 
dTifft k Gregory! 1976; IChincarini k Rood| 1976; prcgory 
\k Thompsorj 1978) show Mayalj had observed foreground 
voidsQ with radi^] ;> 20h 1 Mpc. Joeveer, Einasto, k Tago 



(1978) independently identified voids — holes in their ter- 
minology — in the galaxy distribution in the southern 
galactic hemisphere. 

Other surveys of galaxy types — dwarf, low luminosity, 
irregular, low surface brightness, star-forming, or IRAS — 
that are fo und to respect common voi ds inclu de [Thomp- 
|son] (19 83), iThuan, Gott k Schneider! (1987), |Eder et el 
0989), IBinggeli, Tarcnghi k Sandagej (1990), pothun et 
g] (199jT[ 



Lindcr et al 



(1996), Kuhn, Hopp k Elsasser 
1997), gchombert Pildis fc Eder| (1997), Popcscu, Hop]: 
k Elsasscr| (1997), |Lee et alj (2000), and [El-Ad k Pira H 



(2000). 

The preference of star-forming galaxies for the edges of 
voids certainly agrees with the biased g alaxy formation 
picture, as Salze r has consi stently noted ( (Salzer, Hanson ' 



k Gavazzi 1990; Lee et al. 2000; and references therein) 



Salzer also emphasize s the similar large-scale structures. 



Figure 3 in Lee et al. (2000) clearly illustrates voids de- 



^^The statement in Peebles] (1993) that Kirshner et al 



(1981) 

named these regions voids is wrong; Kirshner et al. clearly state the 
prior discussions. 

2 Hubble's constant is H = lOO/i km s — 1 Mpc -1 . 
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fined by both ordinary and emission-line galaxies. 



Gm clouds also avoid voids. This applies to gas de- 



the exte nsion of the CfA survey (de Lapparent, Geller & 



tected in HI 21-cm emission (iWeinberg et al.| 1991; floff- 



Lu & Salpeter 1992; Bzomoru ct al. 1996; Zwaar 



1997) and gas detected at somewhat lower surface 



man. 
et al 

densities as Lyman-limit or Mgll quasar absorption line 
systems: when the redshift allows the test a galaxy gen- 
erally is observed close to the line of sight and near the 
redshi ft of each absorption line syste m ( Bergeron fc Boissc 
1991; [Bteidel, Dickinson fc Persson| 1994; Lanzctta ct al. 
1995). |Shull, Stocke fc Pcntorj (1996) show that gas clouds 
detected as very low surface density Lyman a absorbers, 
with HI surface density ~ 10 13 cm -2 , avoid dense galaxy 
concentrations. A visual impression is that they also avoid 
the voids, but that awaits a numerical test, perhaps along 
the lines of the nearest neighbor distribution (§3). 

2.2. Void Boundaries 

The strong clustering of galaxies on small scales must 
open low density regions. In a clustering hierarchy model 
that fits the two- through fou r-point correlation functions 
of optically selected galaxies ( [Boneira fc Peebles 1 978) the 
void sizes are compa rable to what is observed ( Soncmi| 



1978; |Vettolani et al.| 1985). The discovery from redshift 
maps is that the void boundaries tend to be rather smooth, 
and defined by galaxies with a broa d ran ge of luminosi- 



ties flChincarini, Rood fc Thompson 1981; de Lapparent 



Geller fc Huchra T986) 



2.3. Void Objects 

Another aspect of the void phenomenon is that isolated 
galaxies do not seem t o be particularl y unusual, apart from 
the gas content. Thus 3zomoru ct al. (1996) conclude that 
"the void galaxies [in BootesJ seem to be unaware of the 
fact that they exist in a huge undcr dense region." In an 
important advance Grogin & Gellei (2000) find evidence 
that the relative velocities of void galaxy pairs at projected 
separation < llbh^ 1 kpc tend to be lower than the relative 
velocities at near cosmic ambient density. Pairs of galaxies 
at separation ~ 100/i _1 kpc typically have other neighbors 
whose massive halos add to the relative velocity, however, 
so the implication for the correlation of halo mass with 
ambient density will require further discussion. 

2.4. Voids and the Morphology-Density Correlation 

The m ix of mo rpholo gical typ es corr elates with ambient 
density ( |Hubble| 1936; |Dresslei] 1980; |Postman fc Gcllcij 
1984). The new aspect of the void phenomenon is that if 
there is a special class of void galaxies the shift to the void 
mix when the ambient density falls below some fraction of 
the cosmic mean has to be close to discontinuous. 



This point seems to have been first made by Kirshner et 
|aT| (1981), who note that if the voids were missing galaxies 
but not mass one might expect low luminosity galaxies are 
more common relative to L ~ L* galaxies in voids. But 
they conclude that if this were so the void galaxies would 
have to be several magnitudes fainter than typical optically 
selected galaxies, because these hypothetical faint galaxies 
are not observed. This agrees with the CfA (C enter for 
Astrophysics) galaxy maps in Figures 2a and 2d of Davis et 
|aL| (1982): the distributions of the high and low luminosity 
galaxies are strikingly similar. The same effect is seen in 



Huchra 1986) 



The morphology-density correlation includes the pref- 
erence of s tar-for ming g alaxies for the edges o f voids 
(iThompsorJ 1983; |Salzerj 1989). But Figure 3 in |Grogin 



fc Gellcrj (2000) shows how subtle is the variation of the 



mix of observed types with ambient density when the den- 
sity is comparable to the cosmic mean. It can be compared 
to the very different mix within voids, if there is a dark 
void population. 

2.5. The Most Extraordinary Objects 

Extremely unusual galaxies are worth special consider- 
ation. The statistics of their environments are insecure, 
by definition, but the discovery of a few cases well inside 
voids would be influential. Two examples of quite unusual 
objects — that happen not to be in voids — may illustrate 
the situation. 

A systematic program of stud y of the remarkabl e blue 
compact ga laxy S BS 0335-052 (|Thuan fc Izotovj 1997; 



2000; iPustiPnik et alT 2000; and references 
shows it has an extended HI cloud with quite 



Vanzi et al 
therein 

low heavy element abundances along the lines of sight to 
star-forming regions, quite young star populations, and a 
large mass-to- light ratio. It has many of the properties one 
might look for in a young galaxy at low redshift, except 
its position. The projected distance to the large spiral 
galax y NGC 1376 is 150 kpc, at close to the same red- 



shift ( [PustiPnik ct alj 2000). Figure 1 shows the position 
of SBS 0335-052 relative to the galaxies in the Optical 
Redshift Survey ( Santiago ct al. 1995)0 with heliocen- 
tric redshifts in the range 4043 ± 300 km s _1 centered on 
SBS 0335-052. This unusual object is not in a dense re- 



3 http:/ /www. astro. princeton. edu~strauss/ors/indox. html 



ii -10 




4 3 
right ascension, hours 



Fig. 1. — Map of Optical Redshift Survey galaxies (filled 
circles) in a slice in redshift space 600 km s _1 deep centered 
on SBS 0335-052 (open box). The field is 23.5h~ 1 Mpc 
wide by 21.3/i _1 Mpc high. 
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gion. More important of the present purpose, it is not in 
a void. The former agrees with the familiar behavior of 
star- forming galaxies. The latter does not naturally agree 
with the biased galaxy formation picture. 

A second example is DDO 154. The larg e ratio of dark 
to luminous mass led |Carignan fc Freemar| (1988) to term 
it a "dark" galaxy. It is not far from the Local Group, 
at supergalactic coordinates SGL= 90°, SGB= +7°, in a 
continuation of the local sheet of galaxies rather than the 
nearby voids. 

A void population would be expected to have a range of 
properties, some with enough stars or gas to be detectable. 
SBS 0335-052 and DDO 154 are strange enough to moti- 
vate the thought that they are unusual members of this 
hypothetical void population, not extremes from known 
populations. But that does not agree with their positions, 
near galaxies. 

These are just two examples. Systematic studies of en- 
vironments of unusual galaxies, particularly early types, 
will be followed with interest. 

3. STATISTICAL MEASURES 

The tendency of less ordinary types of galaxies to avoid 
the voids defined by ordinary spirals has b een probed 



by void p robability funct ions (White 1979; Vogeley et 



aL| 1994; [El-Ad fc Piranj 2000; and references therein) 
two-point correlation functions, and nearest neighbor dis- 
tances. To keep this discussion somewhat limited I discuss 
only the latter two, that have complementary features. 

A particularly useful feature of correlation functions is 
the simple relation between the measurable angular func- 
tion and the wanted spatial function. There is no simple 
relation between nearest neighbor distances in real space 
and the measurable nearest distances in projected angular 
distributions or in redshift space. But as discussed next, 
the interpretation of correlation functions as a probe of 
the void phenomenon is not straightforward. The rest of 
this paper accordingly uses nearest neighbor statistics. 

3.1. Two-Point Statistics 



One sees in Thompson's (1983) maps the tendency of 
Markarian galaxies to avoid both dense regions and voids 
in the distribution of ordinary galaxies. The former means 
that on small scales the Markarian correlation function, 
^MMi and the Markarian-optically selected galaxy cross 
correlation function, £,Mg, are significantly less than the 
galaxy-galaxy function, £ gg . If at larger separations £mm , 
^Mg , and ^ gg had similar values, because both galaxy types 
trace the same large-scale structure, the slopes of £,mm and 
^Mg would be shallower than ^ gg . This would not mean 
Markarian galaxies tend to occupy the voids defined by 
ordinary galaxie s, of course . To the contrary the visual 
impression from Thompson 's maps and the evidence from 



Thompson's nearest neighbor statistic is that Markarian 



galaxies respect the voids. 

The correlation fun ctions for early- and late-type galax- 
ies in Hermit et al. (1996 Figs. 10 and 11) behave as 
just described: at small separations later types have dis- 
tinctly smaller correlation functions, while at separation 
hr ~ 10 Mpc the functions have similar values. This 
certainly demonstrates the small-scale morphology-density 
correlation. One sees from the example of Markarian 
galaxies that it need not imply later-type galaxies tend 



to occupy the voids defined by earlier types. 

If at hr ^> 10 Mpc the correlation function for an un- 
usual type of galaxy were unusually large it could be a 
signature of void galaxies. But because one can trace con- 
nections among occupied regions over quite large distances 
one would also have to consider the way the distributions 
of objects within occupied regions contribute to the abun- 
dances of galaxy pairs at large separations. 

3.2. Nearest Neighbor Statistics 

This measure also is affected by the morphology-density 
correlation, but it may reveal void galaxies through a tail 
in the distribution of nearest neighbor distances. 

The test uses two types of objects. The o-type are ref- 
erence ordinary galaxies, or the proposed equivalent in a 
simulation. The t-type are unusual test objects that may 
have tended to form in the unusual conditions in voids 
between concentrations of o-types. In analyses of obser- 
vations the i-types may be galaxies — dwarf, irregular, 
compact, or low surface brightness — or gas clouds. In 
analyses of simulations the i-types would be mass con- 
centrations that are not expected to develop into ordinary 
galaxies but seem to be capable of forming observable con- 
centrations of stars or gas. 

The distance from a t-type to the nearest o-type object 
is D to , and the distance from an o-type to the nearest 
neighboring o-type is D ao . The probability distribution 
of D to depends on the number density of the reference o- 
type objects, but we have a control from the distribution 
of D oa . If the two types were randomly selected from the 
same underlying population, and the two selection prob- 
abilities differed by a constant factor, the distributions of 
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Fig. 2. — Map of Optical Redshift Survey galaxies in a 
slice in redshift space. The normal to the slice points in 
the direction of the Virgo Cluster, from the opposite side 
of the Milky Way. Filled squares are elliptical galaxies, 
open squares SOs, crosses Sa to Sb-c, plus signs later spiral 
types, and filled circles dwarfs and irregulars. 
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Fig. 3. — Map of ORS galaxies in a slice in redshift space 
on the same side of the Milky Way as the Virgo Cluster. 



D to and D oa would be the same within the noise, even 
if the selection functions varied with position. If the t- 
types tended to be outside the concentrations of o-types it 
would produce a tail of large values in the distribution of 
D to relative to D 00 . 

All applications of this probe have used nearest neighbor 
distances in redshift space, which can be much larger than 
in real space (§3.3). The effect can be analyzed by using as 
separate variables the transverse and redshift separations, 
but this will be left for future work. 



Early examples of this probe are in Boneira & Peebles 
(1977), who used it to test and argue against the idea 
that there is a spatially homogeneous population of field 
galaxies, and Thompson] (1983), who introduced it to the 
analysis of voids 

Nearest neighbor statistics have been appli ed to a con- 



sider able variety of candidate v oid o bjects (Eder et el 
1989: ISalzer, Hanson, & Gavazzil 1990: Rothun ct all 1993 



Pustil'nik et al. 1995; |Lindcr et "aL 1996; Lee et al 
2000). In some cases the nearest reference ordinary galaxy 
is about equally close to test and ordinary galaxies. In oth- 
ers the test sample have the more distant nearest neigh- 
bors. The difference is not very large, however. Thus one 
finds general agreement that all these classes of objects 
avoid the central parts of voids, but mixed opinions on 
consistency with the picture of biased galaxy formation. 
Two new applications of the statistic may help clarify the 
situation. 

3.3. New Examples of the Nearest Neighbor Statistic 

Redshift samples have improved, so an update of the 
nearest neighbor statistic is worthwhile. The example in 
§3.3.1 uses dwarfs plus irregulars as test galaxies, and in 

4 This study was used in §3.1 to illustrate the prnhlrrnat ir interprp- 
tation of correlation functions as a measure of voids. But Thompson) 
(1983) used nearest neighbor statistics, not correlation functions. 



Fig. 4. — Distributions of distances D in redshift space of 
the nearest spiral neighbors of spirals (broken lines) and 
of dwarfs and irregulars (solid lines). The upper dotted 
lines show the effect of randomly shifting one in three of 
the dwarfs and irregulars, the lower dotted lines the effect 
of randomly shifting one in ten. The redshift ranges are 
200 < cz < 500 km s _1 in the bottom plot, 500 < cz < 750 
km s _1 in the middle, and 750 < cz < 1000 km s" 1 in the 
top. 



§3.3.2 low surface brightness galaxies. 

Both exam ples use the spi rals in the Optical Redshift 
Survey (ORS; Santiago et al. 1995) as reference ordinary 
galaxies. At low redshift the ORS broadly samples the 
luminosity function, and it has useful information on mor- 
phology. The authors caution that completeness as a func- 
tion of apparent magnitude and surface brightness varies 
across the sky; the selection may be particularly inhomoge- 
neous for dwarfs and irregulars. The latter need not affect 
the relative distributions of Dto and D oa if the ordinary 
galaxies are homogeneously sampled. And inhomogeneous 
sampling would not seem likely to mask the signature of 
void galaxies. 



3.3.1. Redshift Maps 

Figures 2 and 3 show maps of the ORS galaxies in 
nearby narrow slices in redshift space (corrected to the 
Local Group by adding to the heliocentric redshift 300 
km s _1 toward I — 90° 6 = 0). The normals point to 
supergalactic coordinates SGB — 0, SGL = 90°, roughly 
toward the Virgo Cluster. The slice in Figure 2 is on the 
opposite side of the Milky Way from the Virgo Cluster, the 
slice in Figure 3 on the same side. Galaxies closer than 6° 
from the Virgo Cluster are removed. The slice toward 
Virgo is thinner, because this side is more crowded, and 
it is seriously distorted by virgocentric flow, but early and 
late types may be similarly affected by large-scale stream- 
ing. Most parts of the slices are at high galactic latitude, 
but the upper right corner in Figure 2 and the lower left 
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Table 1 
Nearest Neighbor Parameters 



Test Sample 


rz ■ a 


rz a 


N 


N t 






H D h a 


H (D 00 ) a 


(D t0 )/(D 00 ) 


(D ro )/(D 00 ) 


Dw/Irr 


200 


500 


20 


63 


0.08 


0.26 


11 


110 ± 17 


0.83 ±0.14 


1.43 ±0.23 


Dw/Irr 


500 


750 


51 


66 


0.08 


0.11 


11 


107 ± 16 


0.69 ±0.12 


1.99 ±0.32 


Dw/Irr 


750 


1000 


100 


109 


0.08 


0.09 


11 


78 ±5 


1.10 ±0.12 


3.09 ±0.26 


LSBa 


200 


3000 


153 


27 


0.02 


0.003 


38 


133 ± 10 


0.92 ±0.14 


2.51 ±0.30 


LSBb 


200 


3000 


70 


43 


0.02 


0.01 


32 


115 ± 10 


0.88 ±0.12 


3.19 ±0.36 


LSBa 


3000 


6000 


160 


43 


0.003 


0.0008 


180 


284 ± 17 


1.10 ±0.12 


1.68 ±0.15 


LSBb 


3000 


6000 


80 


60 


0.004 


0.003 


140 


270 ± 20 


1.03 ±0.10 


1.81 ±0.17 


LSBa 


6000 


9000 


106 


42 


0.0007 


0.0003 


560 


490 ± 40 


1.10 ±0.12 


1.39 ±0.17 



a Unit km s" 1 6 Unit h 3 Mpc" 3 



corner of Figure 3 dip to |6| < 30°, where the ORS may 
be significantly less complete. 

The filled squares are elliptical galaxies (Burstein nu- 
merical morph ological type B NMT< 12, de Vaucouleurs 



type T < -5; |Willick ct al 1997, Table 7). The open 



squares are SOs (cut at BNMT< 112, T < 0), the crosses 
Sa through Sb-c (cut at BNMT< 152, T < 4), the plus 
signs later type spirals (cut at BNMT< 182, T < 7), and 
the filled circles all the dwarfs, irregulars, and other cate- 
gories. 

It is difficult to make out much difference in clustering 
properties of the different morphological types in Figures 2 
and 3. It is easy to find nearly empty regions. The former 
is at least in part due to the absence of rich clusters. The 
latter may also be affected by the very limited samples, it 
agrees with the many other examples reviewed in §2 and 
§3-2. 

3.3.2. Dwarfs and Irregulars 

Results of the nearest neighbor test applied to the ORS 
Dw/Irr galaxies (filled circles in Figs. 2 and 3) and spiral 
galaxies (crosses and plus signs) are shown in Figure 4 and 
Table 1. 

Pairs of ORS galaxies closer than 2 arc min and at low 
redshift tend to have the same or quite similar redshifts, 
a sign they may be the same galaxy, so one is discarded. 
I cut the galaxies at galactic latitude \b\ < 30° and within 
6° of the Virgo cluster (a = 186.6°, S = 13.2°). 

The subdivision in the three bins in redshift reduces the 
effect of variable completeness as a function of redshift, 
and it offers a test of reproducibility. 

For each galaxy, o (spiral) or t (Dw/Irr), in each redshift 
bin, Figure 4 shows the distribution of distances in redshift 
space (relative to the Local Group) to the nearest o-type at 
any redshift greater than 200 km s _1 and within the cuts 
in angular position. The plots are normalized cumulative 
distributions. 

In the two lowest redshift bins the Dw/Irr galaxies tend 
to have the closer spiral neighbor; in the highest redshift 
bin the spirals tend to have the closer spiral neighbor. The 
difference may be related to the variation in the range 
of absolute magnitudes sampled in each redshift bin (re- 
flected in the variation in the ratio of numbers N t and N D 
of Dw/Irr and spirals). Or it may only be noise, for the 
differences in the distributions are small. 



Table 1 shows the means (D 00 ) of the distances of 
the nearest spiral neighbors of spirals, and the ratios 
(D to ) I (D 00 ) of the means of nearest neighbor distances 
of Dw/Irr and of spiral galaxies. The standard devia- 
tions assume each distance is statistically independent. (If 
two o-types are each other's nearest neighbor the distance 
counts once.) The differences of the distributions in Fig- 
ure 4 may appear more significant than the differences of 
the ratios (D to )/ (D 00 ) from unity, because the cumulative 
distributions over-emphasize the noise. 

We can estimate the nearest neighbor distance among 
spirals in real space from the observation that the small- 
scale distribution of optically selected galaxies approx- 
imates a scale- invariant clustering hierarchy, or fractal, 
with two-point correlation function £ = (r G /r) 7 . The dis- 
tance rh at which a spiral has on average one spiral neigh- 
bor satisfies 



Cd 3 r _ AlTn ° r l 

.! 3-7 



= 1, 



(1) 



if the mean density n is large enough that rh <C r Q . Since 
the small-scale distribution is close to scale-invariant the 
mean Dh of the nearest neighbor distance differs from rh 
by a fixed factor. Since the factor is not likely to be greatly 
different from unity a useful estimate of the mean real 
distance from a spiral to the nearest spiral is 



3-7 
47rn rJ 



1/(3-7) 



(2) 



The table lists H Q D h for hr a = 5 Mpc and 7 = 1.77. 

The small-scale relative peculiar velocity in the field is 
an order of magnitude larger than H a Dh in the ORS at 
cz < 3000 km s -1 . This means the ratio {D to ) / (D QO ) 
is determined by peculiar velocities; it does not tell us 
whether Dw/Irr or spiral galaxies tend to have the closer 
spiral neighbors in position space. 

The D to could be reduced by accidental cancellation of 
differences in the distributions of relative velocities and 
positions, but that seems unlikely. Thus I conclude 

(1) the distributions of peculiar velocities of Dw/Irr and 
spiral galaxies relative to nearby spirals are quite similar, 
and 

(2) the mean distance in position space from a Dw/Irr 
to the nearest spiral is less than about l/i -1 Mpc. 
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3.3.3. Low Surface Brightness Galaxies 




1000 

D, km/s 



1500 2000 



Fig. 5. — Distributions of distances of the n earest spiral 

" (1992) LSB 



neighbors of spirals and of Schombert et al. 
galaxies. The dotted lines show the effect of randomly 
shifting one in three of the LSB galaxies. The redshift 
ranges are 200 < cz < 3000 km s _1 in the bottom plot, 
3000 <cz < 6000 km s" 1 in the middle, and 6000 < cz < 
9000 km s" 1 in the top. 



The second result depends on the number density of spi- 
rals, here n ~ O.lh 3 Mpc -3 . But, since most spirals avoid 
the voids, if n a were decreased by cutting the sample at 
larger luminosity or circular velocity, the Dw/Irr's would 
still tend to be within 1ft, -1 kpc of the concentrations of 
the sample spiral galaxies. 

The last column in the table shows the result of moving 
each of the Dw/Irr galaxies to a randomly chosen position, 
with constant probability per unit solid angle within the 
cuts in galactic latitude and distance from the Virgo Clus- 
ter, and at uniformly probability distribution in z 3 within 
the redshift bin for the galaxy. The increase in nearest 
neighbor distances is significant, though not all that large. 

The separate test illustrated by the upper dotted lines in 
Figure 4 shows the result of shifting to a random position 
one in three of the Dw/Irr (selecting every third object 
in the catalog list for a new position), leaving the rest at 
their catalog positions. The lower dotted lines show the 
effect of randomly repositioning one in ten. At one in 
three the resulting tail in the nearest neighbor distances 
is marginal in the shallowest sample, pronounced in the 
two deeper ones. It follows that at the catalog positions 
fewer than one in three of the Dw/Irr's could be members 
of a homogeneously distributed population, for otherwise 
the distributions of D to for the other galaxies could not be 
tight enough to match the observed distribution of all D to - 
When one in ten is randomly shifted the tail is pronounced 
only in the deepest sample. Thus I conclude 

(3) fewer than about one in ten of the Dw/Irr is in a sta- 
tistically homogeneous randomly distributed population. 



The Schombert et al. (1992; here LSBa) and Impey et 



(1996; here LSBb) low surface brightness galaxies are 
analyzed separately to check reproducibility. The refer- 
ence ordinary galaxies are the ORS spirals used in §3.3.2. 

Figure 5 shows the nearest neighbor distributions for the 
LSBa and ORS spiral galaxi es at declination 10 < 6 < 25°, 
which includes most of the Schombert et al.| (1992) low 
and very low surface brightness galaxies, at galactic lat- 
itude \b\ > 30°, and in the indicated ranges of redshift 
corrected to the Local Group. To be counted the near- 
est spiral neighbor must be more than 2 arc min away, 
to reduce the chance the same object is in both catalogs, 
and at cz > 200 km s -1 , < S < 35°, and \b\ > 30°. 
The cuts for LSBb (Fig. 6) are the same except for the 
declination: the galaxies whose neighbors are counted are 
at —1.5° < S < 3.5°, and the spiral neighbors are at 
-11.5° < S < 13.5°. The LSBb galax ies are the la rge an- 
gular size objects with HI redshifts in [mpey et al. (1996). 
(The optical redshifts are not used because their uncer- 
tainties would significantly add to the nearest neighbor 
distances in redshift space.) 

At the mean number density n of spirals in the nearest 
redshift bin, 200 < cz < 3000 km s _1 , the characteristic 
physical distance to the nearest neighbor in the clustering 
hierarchy is H Dh ~ 35 km s _1 . As for the Dw/Irr case, 
this is small compared to typical relative peculiar veloc- 
ities. The quite similar distributions of distances in red- 
shift space to the nearest spiral neighbors of these LSBs, 
Dw/Irr's, and spiral galaxies, imply both LSB and Dw/Irr 
galaxies typically are closer than 1ft -1 Mpc from a spi- 
ral and all three galaxy types have quite similar relative 
peculiar velocities. 

In the deepest redshift bin (the bottom line in Table 1) 



T3 

a 

o 



a. 




0.01 



1000 



D, km/s 



Fig. 



The same as Fig. 5 for the Impey et al 



LSB galaxies at redshift ranges 200 < cz < 3000 
in the bottom plot and 3000 < cz < 6000 km s _1 
top. 



(1996) 
km s _1 
in the 
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the mean nearest neighbor distance (D ro ) for randomly 
placed LSBa's is not significantly larger than the mean 
(Dto) at the catalog positions. This can be understood as 
follows. In a homogeneous random Poisson process with 
the mean number density n Q = 0.0007ft 3 Mpc -3 of this 
sample of spirals the mean distance to the nearest neighbor 
is 

r(1/3) - 6ft- 1 Mpc. (3) 



D 



Pois 



(36 7 rn) 1 /3 



This is comparable to (D 00 ) and to the clustering length 
r D . This means the sparse sampling has suppressed the 
sensitivity of the nearest neighbor statistic to the galaxy 
distribution. Figure 6 for LSBb accordingly shows only 
the two lower redshift bins. 

At the number density of spirals in the bin 3000 < cz < 
6000 km s _1 the physical distance Dh ~ 150 km s^ 1 in the 
clustering hierarchy is not much less than the relative ve- 
locity dispersion. The physical separations thus might be 
expected to make a significant contribution to the nearest 
neighbor distances in redshift space. Consistent with this, 
(D 00 ) is not much larger than Dh- 

The dotted lines in Figures 5 and 6 show the effect of 
randomly moving one in three of the LSB galaxies, leaving 
the rest at their catalog positions. The effect is marginal 
in some cases, but reproducible enough to show that fewer 
than one in three of the LSBs could be in a statistically 
homogeneously distributed population. Since most of this 
hypothetical homogeneous population would be in voids I 
conclude that 

(4) the number density of observed LSB galaxies satisfies 



71LSB < 0.001ft 3 Mpc" 



(4) 



— BoLhuil el al. (1993) find that the nearest spiral neigh- 
bor of an LSB is on average 1.7 times further than that of 
a spiral galaxy. The difference from the results presented 
here is a cautionary example of the sensitivity to samples 
and methods of analysis. 

Bearing in mind this example, but considering also the 
consistency of the distributions in Figures 5 and 6 from 
two independent low surface brightness samples and two 
redshift bins, I conclude that 

(5) the mean nearest neighbor distances in redshift space 
for spirals and LSBs are not likely to differ by more than 
about 30 percent. 

Since this analysis uses low surface brightness galaxies 
detected in HI the result is in line with the HI surveys in 
emission and absorption that show that gas clouds very 
distinctly prefer to be near galaxies (§2.3). But the quan- 
titative constraints may be useful. 

3.4. Summary of the Nearest Neighbor Measure 

The candidate ordinary optically selected galaxies in a 
simulation of galaxy formation ought to define realistic 
voids. If the mean nearest neighbor distance among the 
candidate galaxies scales according to equation (0), 



hD h 



0.6 



0.01ft 3 

n Mpc 3 



0.8 



Mpc, 



(5) 



when the cut in luminosity function or circular velocity 
produces mean number density n 01 the distribution may 



be expected to approximate the observed scale-invariant 
clustering hierarchy. In this case the voids likely are sat- 
isfactorily large. 

Candidate void galaxies in a simulation of galaxy forma- 
tion are the mass concentrations that are not good candi- 
dates for ordinary galaxies — perhaps their circular veloc- 
ities are too small — but seem to be capable of developing 
into observable gas clouds or galaxies of stars. They might 
be candidates for Magellanic-type irregulars, with circu- 
lar velocities greater than about v c ~ 20 km s _1 , large 
enoug h to r esist substantia l loss o f photoion ized plasma 



(Rees 1986; Babul k Rees 1992; Efstathiou 1992). Or 



they might be candidate low surface brightness or com- 
pact galaxies. 

In the examples in §3.3 the mean physical distance from 
a Dw/Irr or LSB galaxy to the nearest spiral satisfies 

(Ao)physical £ lft _1 MpC, (6) 

when the mean number density of spirals is 

n a £ 0.01ft 3 Mpc~ 3 . (7) 

The samples avoid rich clusters, that would increase the 
mean distances in redshift space and reduce them in real 
space. In redshift space equation (||) is replaced with an 
equality to about 10 percent. 

At smaller n the examples in §3.3 indicate 



(-^£o)physical ^ 2(Z} GO )physical ^ 2_Dfr 



(8) 



The factor of two takes account of applications of the near- 
est neighbor test that find (D to ) > (D 00 ). In the examples 
presented here (D to ) and (D 00 ) are equal to ~ 20 percent. 

These numbers can be considered a quantitative defini- 
tion of the void phenomenon. One must add the cosmic 
number density, nt, because a model may satisfy equa- 
tion (|J) with a high density of void galaxies and an unac- 
ceptably high density of clustered objects. Equation (^) 
gives a limit on the density of candidate void LSB galax- 
ies; the corresponding limit on void dwarf plus irregular 
galaxies is 

^Dw/Irr £ 0.01ft 3 Mpc^ 3 . (9) 

More observational checks of these numbers would be 
of interest. The application to realizations of models for 
galaxy formation might interesting too. 

4. THEORETICAL SITUATION 

The conventional theoretical interpretation of voids was 
driven by the elegance of the Einstein-de Sitter model 
(with fi m = 1 in matter capable of clustering, and neg- 
ligibly small space curvature and cosmological constant), 
for this cosmological model requires most of the mass to 
be in the voids. It also is motivated by the CDM model 
for structure formation, which naturally produces a biased 
distribution of galaxies relative to mass. The considera- 
tions prior to the paradigm shift to f2 m = 0.25 ± 0.1 are 
worth reviewing as a guide to the present situation. 

4.1. Voids in an Einstein-de Sitter Universe 

The small relat ive velocity dispersion in the CfA sample 
( Davis k Peebles! 1983) shows that if fl m = 1 then most of 
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the mass h as to be in the voids Davis et al. (1985; follow- 
ing Kaiser 1984 and Bardccn 1986) show this could occur 
in a natural way in the CDM model: if ordinary galaxies 
formed preferentially in high density regions they would be 
strongly clustered, leaving most of the mass in the voids. 
The void mass would be clumpy and thus might be ex- 
pected to produce bound objects, some observable though 
different fro m ordin ary galaxies. This seems intuitively 
reasonable ( Peebles 1986, 1989). I t is mad e quantitative 



in com putations by Einasto et al. (198 4), | Dckel fc Silk 



Villumsenl (1992) and |Hoffmaii, Sil~ 



(1986), |Braine rd 
Wyse|(1992)J3 

Opinions on whether the distributions of galaxy types 
argue for or against the CDM model are mixe d. Argu 



mcnts in favor cite morphological segregation (Dekel & 
SilkfT|986), the smaller two-p oint correlation function for 
classe s of unusual galaxies ( Salzcr, Hanson, fc Gavazz] 



contain little mass. This is the natural interpretation of 
the void phenomenon. But this consideration did not play 
a significant role in the change of the most favored model 
from Einstcin-de Sitter to low Q m . I suspect one reason is 
the prediction by CDM simulations that the voids contain 
significant mass even when fi m is small. 

A visual impression of numerical simulations of the low 
density ACDM model (with a cosmological constant to 
make flat space sections) is that of classical biasing: the 
larger dark mass halos cluster more strongly than the 
mass, and the less massive halos spread into the voids de- 
fined by the larger halos. It is no criticism of these studies, 
which report the behavior of the model, to note that one 
is not reminded of the void phenomenon illustrated in Fig- 
ur es 2 and 3. The sit uation is clearly presented in Figure 5 
~ (1999). 



of Kauffmann et al 



1994), an d low surface 



1990; [Mo, McGaugh fc Bothun 

brightness galaxies ( |Hoffman, Silk fc Wyse 1992). Section 
3 presents reasons for treating these lines of evidence with 
some caution. Others conclude that the void pheno menon 



The remedy may be a proper understanding of how 
the form ation of o bserva ble ob jects may be suppressed 
in voids (Ostriker 1993). Rees (1985) explores the pos- 



challenges biased galaxy fo rmat ion ( Einasto et al. 1984 



1989; 
al 



Binggeli, Tarenghi & Sane 



age 



19)4 



1990; and Einasto ei 



sible role of relativistic particles or ionizing radiation from 
the first generation of galaxies in suppressing subsequent 
galaxy formation, that otherwise would tend to occur in 



Thuau, Gott fc Schneider] 1987; Eder et el. 19 89; Peebles voids. Cen & Ostriker (2000) analyze realizations based 



Ostrikor 1 (1993) gives a balanced summary of the gen- 
eral opinion in the community shortly before the paradigm 
shift to low Q m ; "Nominally one expects, in the CDM 
model, that the voids will be populated to a degree larger 
than is observed ... but in the absence of agreed-upon 
theories of galaxy formation, it is difficult to quantify this 
apparent disagreement." 

Ostriker's assessment is accurate, but we do have guid- 
ance on what might have happened from what is observed 
( Peebles] 1989). Consider the Magellanic-type irregulars 
on the outskirts of the Local Group — IC 1613, Sextans A 
and B, WLM, IC 5152, and NGC 3109 — at distances 



on a numerical prescription for how galaxies form or are 
suppressed. They find the prescription produces obser- 
vationally acceptable realizations. An application of the 
nearest neighbor statistic in §3.4 would be interesting. An 
explanation of the Im phenomenon (§4.1) would be edify- 
ing. 

The remedy may be a n adjustment of the CDM model. 
Bode, Ostriker & Turok (2000) show that warm dark mat- 



between 0.7 and 1.7 Mpc ( |Tullyj 1988; Table 1 in |Peebles 
et al.| 2000). They have small peculiar velocities relative 



to the Local Group. They are not near a large galaxy, so 
they are not likely to have been spawned by tidal tails or 
other nonlinear process in the large galaxies. Since they 
are at ambient densities close to the cosmic mean their 
first substantial star populations would have formed un- 
der conditions not greatly different from the voids at the 
same epoch. Why are such galaxies so rare in the voids? 



ter can produce regions that are quite devoid of gravita- 
tional seeds for structure formation. The effect is strik- 
ing, but the scale unfavorable. If in their Figure 11 the 
large circles represented ordinary optically selected galax- 
ies, which would make about the observed number density, 
and the small circles were dwarfs or irregulars, the model 
would not seem to agree with the nearest neighbor distri- 
bution in Figure 4, or the maps in Figures 2 and 3. If 
the scale were enlarged so all circles represented ordinary 
L ~ L* galaxies, the voids would be well represented, but 
the Lyman-a forest at z = 3 would be a problem. 

5. SUMMARY REMARKS 
5.1. The Observational and Theoretical Situations 



If fL 



4.2. Voids in a Low Density Universe 
= 0.25 ±0.1 the observations are consistent with 



1. Some galaxy types prefer dense regi ons (Hubble 
1936), others the edges of voids (Thompson 1983). This 



the assumptio n that ordinar y optically selected galaxies 
trace the mass (Bahcall et al. 2000 and references therein). 
If galaxies are good mass tracers we can assume the voids 

5 The relative velocity dispersion in the CfA sample is biased low 
by the i i nder-represent ati on of rich clusters with large velocity dis- 
persions (fVIarzke et al. 1995). But the mass within the Abell radii of 
the Abell clusters is only about one percent of the critical Einstein- 
de Sitter value, and the low relative velocities of galaxies outside the 
clusters indicates the mass in the less dense galaxy concentrations 
sampled by CfA is significantly less than critical too 

6 The COSmic String picture with hflt r\xr\r maftnr f^r-hm-rpr oi 



1989), and the explosion variant (Ostriker, Thompson & Wit 



tern 1986), might produce voids that never were substantially dis 
turbed, though that does not agree with the recent picture of a near 
space-filling Lyman-a forest at redshift z ~ 3. 



well-established morphology-density correlation ( Dressier] 
1980) seems to arise in a natural way from hierarchical 
gravitational structure formation, as in the CDM model 
( [White et al] 1987). 

2. All known galaxy types and most gas clouds are 
scarce outside the concentrations of ordinary galaxies. 
This void phenomenon is discontinuous from the morphology- 
density correlation (§2). It appears to be a more challeng- 
ing test of ideas on how galaxies formed (§4.2). 

3. Two-point correlation functions are sensitive probes 
of the morphology-density correlation but, I argue in §3.1, 
are not readily interpreted measures of the void phe- 
nomenon. 

4. The advantage of the nearest-neighbor statistic is a 
reasonably simple interpretation as a constraint on void 
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galaxies. This is widely recognized in the observational 
community (§3.2). The summary in §3.4 may be useful 
for testing simulations of galaxy formation. 

5. The void phenomenon is observed among an impres- 
sively broad range of objects. The weight of this evidence 
naturally is from more readily observable gas-rich and star- 
forming objects, however. It m ay be significant that the 
considerations of Dekel fc SU3 (1986) argue for early-type 
void objects. Further observational tests for such objects 
are particularly desirable. 

6. The challenge to the ACDM model might be resolved 
by a demonstration that the formation of observable void 
objects really can be adequately suppressed. This ap- 
proach is challenged in turn by the Magellanic-type ir- 
regulars on the outskirts of the Local Group, that seem to 
have formed in ambient conditions not very different from 
the voids (§4.1). 

7. The challenge might be resolved by adjusting the 
model for structure formation. A perhaps desperate idea 
is that the voids have been emptied by th e gravit ational 
growth of holes in the mass distribution (Peebles] 1982). 
The void phenomenon seems striking enough to motivate 
a search for viable initial conditions for this picture. 

5.2. Interpretations of Voids as an Example of the 
Scientific Method 

Th e introductio n of simulations of biased galaxy forma- 
tion ( Davis ct al.| 1985) was not inspired by or even ob- 
viously consistent with the evidence that giant and dwarf 



galax ies have qui te similar distributions (Kirshncr ct al 
1981; Davis et al. 1982). This follows an honorable tradi- 
tion in cosmology and, I suspect, other developin g sciences . 
A strikingly successful example in cosmology is Einstein 's 
(1918) cosmological principle. It did not agree with what 
was then known, but it led us to an aspect of physical 
reality. 

The advances in ideas on structure formation since 
1985 have not been seriously influenced by the void phe- 
nomenon. Again, this is not unusual. Another example is 
the galaxy n-point correlation functions. The near power 
law form of the observed two-point function is a widely 
discussed test. The three-point function is little noted in 
discussions of simulations, despite its importan ce in char 
acteri 



ing the amall acalc galaxy diatribution ( S- 



aoncira 



3" 



Peebles 1077). — Thie complex and eomotimoe weak intor 



pla y of th eory and practice in cosmology is well represented 
by Kuhn's (1962) paradigms, with socially selected theo- 
ries and constraints. The weak interplay can be healthy. 
It allows a concentrated study of a particular subset of 
ideas, that may establish or eliminate them. It reduces 
the chance of distraction by misread evidence. It also al- 
lows distraction by unprofitable ideas, of course, but that 
is remedied when theory and practice mature and on oc- 
casion produce crises that drive paradigm shifts to better 
approximations to reality. 

The CDM model is maturing, most dramatically in its 
success in relating the power spectrum of the thermal back- 
ground radiation temperature to observationally accept- 
able cosmological parameters (eg. Hu et al. 2000). This 
shows the CDM model likely is a good approximation to 
how structure started forming on the length scales probed 
by the measurements. 

The apparent inconsistency between the theory and ob- 



servations of void is striking enough to be classified as a 
crisis for the CDM model. It may be resolved within the 
model, through a demonstration of an acceptable theory 
of galaxy formation. Or it may drive an adjustment of the 
model. 

I have benefitted from discussions with Jerry Ostriker, 
Stacy McGaugh, Michael Strauss, Laird Thompson, Trinh 
Thuan, Brent Tully, and Neil Turok. This work was sup- 
ported in part by the National Science Foundation. 
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